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Abstract 

The possibility of generating multiple coherent modes in trapped Bose gases is ad- 
vanced. This requires the usage of several driving fields whose frequencies are tuned 
close to the corresponding transition frequencies. A general criterion is derived ex- 
plaining when the driving fields, even being in perfect resonance, cannot generate the 
topological coherent modes. This criterion is termed the theorem of shape conserva- 
tion. Bose-Einstein condensates with generated coherent modes display a number of 
interesting effects, such as: interference fringes, interference current, mode locking, dy- 
namic transition, critical phenomena, chaotic motion, harmonic generation, parametric 
conversion, atomic squeezing, and entanglement production. Approximate solutions, 
based on the averaging techniques, are found to be in good agreement with direct 
numerical calculations for the Gross-Pitaevskii equation. 
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1 Introduction 



Dilute Bose gases at low temperatures are described by the Gross-Pitaevskii equation (see 
reviews [1-6]). The ground-state solution of the stationary Gross-Pitaevskii equation corre- 
sponds to the equilibrium Bose-Einstein condensate. In the presence of a trapping potential, 
this equation possesses a discrete spectrum and, respectively, a discrete set of wave functions 
describing the topological coherent modes, which represent nonground-statc condensates [7]. 
These modes are called topological since the related wave functions have different spatial 
shapes, with different number of zeros. And the modes are termed coherent because the 
related wave functions correspond to coherent states [8]. Vortices are a particular case of 
such modes [9,10]. 

In practice, the topological coherent modes can be generated by subjecting the system 
of trapped atoms to the action of an alternating external field with a frequency tuned to the 
resonance with a transition frequency between the ground state and an excited collective 
level [7,9-12]. The properties of these modes have been studied in several pubhcations [7- 
26]. Instead of modulating with a resonant field the trapping potential, one can periodically 
vary the atomic scattering length by invoking the Feshbach resonance [27-29] . Note that the 
coherent modes, being the solutions to the nonlinear Gross-Pitaevskii equation, in some cases 
can be considered as analytical continuations, under increasing nonlinearity, of the related 
linear modes; but there exist also the irreducible modes having no linear counterparts. The 
latter happens for sufficiently strong nonlinearity due to atomic interactions, when atoms are 
in double- or multiwell potentials [19,22] or inside an optical lattice [30-33]. Thermodynamics 
of the coherent modes has also been considered [34]. 

One should not confuse the topological coherent modes, which are the solutions to the 
nonlinear Gross-Pitaevskii equation, with the elementary excitations treated by the linear 
Bogolubov - De Gennes equations. In principle, one could consider resonant transitions 
between the modes of elementary excitations [35], which, however, is a radically different 
problem. 

In previous publications, the case of a single generated coherent mode has been studied, 
which requires the action of one external resonant field. Here, we generalize the consideration 
to the multiple mode generation, which can be done by involving several alternating fields. 
It is particularly important that we have accomplished all investigations in two ways, by 
employing approximate analytical methods, based on the averaging techniques, and also by 
solving numerically the Gross-Pitaevskii equation, both ways being in good agreement with 
each other. 



Topological coherent modes are the solutions to the stationary Gross-Pitaevskii equation 



2 Multiple Mode Generation 



(1) 



with the nonlinear Hamiltonian 



V'^ + U{r) + NA,\(p 



(2) 



2mo 
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containing a trapping potential U{r) and the nonlinear term due to the binary contact 
interactions 

A, = Anh^ ^ ; 

mo 

N being the number of particles in the trap; mo, particle mass; a^, scattering length. The 
functions v'n(r) are normalized to unity, (<^„, (fn) — 1- In general, and are not 
compulsorily orthogonal. 

Temporal evolution is described by the time-dependent Gross-Pitaevskii equation 

th cp{r,t) ^ {h[<p] + v) <p{r,t) , (3) 

with an external potential V — V{r, t), which is assumed to have the multifrequency form 

^ = ^ E [BA^)e'^'' + B;{v)e-^^^'] , (4) 

j 

in which j = 1,2, . . .. The frequencies of the alternating potential (4) are chosen so that 
each of them is tuned close to one of the transition frequencies 

<^mn = ^ (Em — En) (5) 

between two coherent modes, which implies the validity of the resonance conditions 

< 1 , ^mn = l^j - l^mn ■ (6) 

Transitions between different modes can be produced by the action of the multiresonant 
potential (4) and also owing to interatomic interactions. Consequently, there are two types 
of transition amplitudes, defined as the matrix elements over the related coherent modes. 
The transition amplitudes 

Oimn = N ^ {\ipn\^ , 2\ipn\'^ - \ip„,\'^j (7) 

correspond to interatomic interactions, while the transition amplitudes 

Pmn = ^ ((Pm,Bj(fn) , (8) 

with Bj = Bj(r), are characterized by the related resonant parts of the potential (4). Good 
resonance can be supported provided the transition amplitudes (7) and (8) are much smaller 
than the related transition frequencies. 



Or 



«1, 



< 1 . (9) 



The first of these inequalities limits the number of atoms in the trap [12]. This limiting 
number of atoms is of the order of the critical number of atoms with attractive interactions 
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allowing for the stability of the trapped atomic gas [7,12,36]. Hence, the generation of coher- 
ent modes is admissible for atomic systems with repulsive as well as attractive interactions. 
The solution to Eq. (3) can be presented as the expansion 

ip{r, t) = Y^ Cn(t)<^n(r) exp ^ Ent^ (10) 

over the coherent modes. The coefficient functions c„(i) are assumed to be slowly varying 
in time as compared to the exponentials, so that 



dCr 



dt 



« 1 . (11) 



Inequalities (8) and (11) are mutually self-consistent, since if a^n and f3rnn ^ 0, then 
Cnit) ~^ const and dcn/dt — > 0. Substituting expansion (10) into the evolution equation 
(3) and employing the averaging technique [7,12] results in the system of equations for the 
functions Cn{t). The values |cn(t)p describe fractional mode populations. For instance, in 
the case of three modes coupled by two resonant fields, we find 

I ^ \Oii2\C2\ +Q;i3|C3| jCi + fi, 

dC2 

'dt 

dC3 

'dt 

where the terms fj depend on the kind of the mode generation involved. Thus, for the 
cascade generation. 



• "^2 / I |2 , I |2\ , 

I ^ [a2l\Cl\ +Q;23|C3| )C2 + f2, 

^ -IT = («3l|ci|^ + Q!32|C2|^) C3 /3 , (12) 



fi = l /3i2C2e^^-* , f2 = l /3r2Cie-*^-* + ^ /523C3e^^-* 



f3 = l /323C2e-'^-* ; (13) 



for the y-type generation. 



fs^l /9i*3Cie-'^-* ; (14) 



and for the A-type generation 



h^l Plsc.e-^^''' + \ /?23C2e-^^-* . (15) 
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Equations (12) are to be complimented by initial conditions and by the normalization con- 
dition 

|ci|' + |c2p + leap = 1 . 

Note that Eqs. (12) enjoy the global phase symmetry Cj — * Cj exp(ix), with x ^ [0, 27r]. This 
could allow us to take one of the initial conditions for cj (0) as a real quantity. However, one 
cannot make all three functions Cj{t) real. In general, they are complex- valued functions. 

Let us emphasize that all generated modes coexist in the same trap and are not spatially 
separated. This distinguishes the considered situation, both physically as well as mathe- 
matically, from the cases of spatially separated condensates, as in multiwell potentials [37] 
or in clouds separated by means of the collective Rayleigh scattering [38-41]. The coherent 
modes, we consider, differ from each other by the shapes of their wave functions and by their 
energies, which are defined by the Gross-Pitaevskii equation (1). 

3 Theorem of Shape Conservation 

Applied alternating potentials not always can generate new topological coherent modes, even 
if the required resonance conditions are perfectly valid. It may occur that the initial atomic 
cloud, being subject to the action of resonant fields, oscillates in space without changing its 
shape, but no higher modes are excited. The criterion for such a behaviour is provided by 
the theorem below. 

Let us consider Bose-condensed atoms in a trapping potential U{r), with a driving field 
V{r,t). The corresponding coherent wave function (p{T,t) satisfies the temporal Gross- 
Pitaevskii equation (3). Since atoms are trapped, then the trapping condition 

V9(r,t)->0 (|r|^oo) (16) 

is always valid for all t > 0. Assume that at the initial time the function ip{r, 0) represents 
a real mode 

<p{v,0)^<po{r)^<pl{r) (17) 

satisfying the stationary equation (1). We shall say that the density of atoms conserves its 
shape when the shape-conservation condition 

|(^(r,t)| = |(^(r-a,0)| (18) 

holds true, where a = a.{t) defines the center-of-mass motion. The following criterion has 
been proven [42]. 

Theorem. The solution to the Gross-Pitaevskii equation (3), obeying conditions (16) 
and (17), satisfies the shape-conservation condition (18) if and only if the trapping potential 
is harmonic, while the driving field is linear in space, that is, 

U{r) = A + Ai • r + J2A^pr''r^ , V{r) = Boit) + B^{t) ■ r , (19) 

where -Bo(^) Bi(t) are arbitrary functions of time. Then the center-of-mass motion is 
given by the equation 

^0 ^ + E(^«/J + Af,^)a^ + B^{t) = . (20) 
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This theorem shows that, even when the functions -Bo(^) and Bi(t) in Eq. (19) correspond 
to alternating fields being in resonance with some transition frequencies, no generation of 
other modes is possible, but the density of atoms will conserve its shape, with the ccnter-of- 
mass motion described by Eq. (20). Hence, in order to effectively generate coherent modes, 
one has to avoid condition (19). 



By applying k alternating fields, with the frequencies cuj {j — 1,2, ... ,k) one can generate 
k topological modes, starting with one given mode, say, having at the initial time all atoms 
in the ground state. Then we will get A; + 1 coupled coherent modes. The mode ampli- 
tudes Cj, defining the fractional populations |cjp, obey k + 1 nonlinear differential equations. 
Because of their nonlinearity, these equations display a rich variety of interesting features. 
We have thoroughly examined the dynamics of Cj for the cases of two and three coupled 
coherent modes, finding all fixed points and accomplishing a complete stability analysis. For 
the corresponding physical cases, we have also solved the Gross-Pitaevskii equation numer- 
ically. Both approaches were found to be in good agreement with each other. Since the 
derivation of the equations for Cj is based on the averaging technique [43], the resulting 
equations are approximate, with the accuracy of their solutions characterized by the values 
of the involved parameters. The errors of such approximate solutions are of the order of 
ma.x{\amn/'-^mn\, |/5mn/'^mn|}- Within this accuracy, the results obtained by solving the evo- 
lution equations for the mode amplitudes Cj coincide with those following from the direct 
numerical solution of the Gross-Pitaevskii equation. Details of these calculations will be 
published in a separate paper. And here, being limited in space, we shall present a brief 
account of the most interesting physical effects we have found. 

(1) Interference Fringes 

The total density of atoms in the trap, p{r,t) = \(p{r,t)\'^, is not equal to the sum of 

the partial densities pj{r,t) = |(/9j(r, t)p where ipj{r,t) = Cj{t)ipj{r), but there exist the 
characteristic interference fringes described by the interference density 



The appearance of such interference patterns is typical of coexisting coherent states [12,20,25,26]. 

(2) Interference Current 

Because of an essential nonuniformity of the atomic density for several coexisting coherent 
modes, there arises the interference current 



which, similarly to the interference density (21), is the difference between the total current 



4 Dynamic Resonant Effects 



Pint{r, t) = p(r, t)-Y^ Pji^^, t) ■ 



(21) 



j 




(22) 



j(r,t) 



h 



Im (p*{r,t)V(p{r,t) 



mo 



6 



and the sum of the partial currents 



ji{r,t) = — Im ip*{r,t)Wipi{r,t) . 
mo 

The interference current is also typical of the case of different coexisting coherent modes 
[12,20,25,26]. Such a current is called sometimes the internal Josephson current or topological 
current. The possibility of the Josephson-type oscillations between two interpenetrating 
populations, not separated by any barrier, was suggested by Leggett [44]. 

(3) Mode Locking 

When the amphtudes of the driving resonant fields are sufficiently small, such that 
iPmn/oimnI <S 1, the fractional mode populations exhibit nonlinear oscillations close to their 
initial values, never crossing the line 1/2, that is, depending on initial conditions, one has 
either 

< |c,r < I (23) 

or 

1 < i^r < 1 . (24) 

The mode locking effect exists for two [7,11,12] as well as for several [42] coupled modes. 
Mathematically, it is the same as self-trapping [45] of atoms in one of the wells of a stationary 
multiwell potential. 

(4) Dynamic Transition 

Increasing the amplitudes of the driving fields up to the values such that \Pmn/(^mn\ ~ 0.5, 
one passes from the mode-locked regime of motion to the mode-unlocked regime, when the 
mode populations oscillate in the whole region between zero and one, 

< |cjf < 1 . (25) 

The change of this dynamic behaviour occurs when the starting point of a trajectory is 
crossed by a saddle separatrix [11,12,25,26]. A similar effect of the change from the self- 
trapped to untrapped regime exists for multiwell potentials [46]. 

(5) Critical Phenomena 

On the manifold of the system parameters Pmn, and A^n, there exist surfaces 

whose crossing results in the dynamic transitions from the mode-locked to mode-unlocked 
regimes of motion. In the vicinity of these separating surfaces, the temporal behaviour 
of the mode populations experiences sharp changes [11,12,18,20,26], because of which the 
separating surface can be termed the critical surface. Moreover, the time-averaged system 
exhibits on a critical surface critical phenomena analogous to those happening in equilibrium 
statistical systems under second-order phase transitions [20]. 

(6) Chaotic Motion 

When the amplitudes of the driving fields become so large that the ratio \l3mn/<^mn\ is 
of order one, the temporal evolution of the fractional mode populations |cjp can go chaotic. 
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This takes place only when the number of the coupled coherent modes is equal or larger 
than three. Two coexisting modes display always a regular behaviour [11,12,26]. The onset 
of chaos for three or more modes is connected with the disappearance of stable fixed points 
[42]. This is similar to the arising chaos in a three- well potential [37]. 

(7) Harmonic Generation 

In order to generate a coherent mode, it is not compulsory to invoke the standard res- 
onance, as in Eq. (6), when the frequency uj of the driving field is close to the required 
transition frequency UJ21 = (£'2 — Ei)/h. But, instead of the resonance condition u = LJ21, 
one can employ one of the harmonic generation conditions 

nu! = UJ21 (n = 1, 2, . . .) . (26) 

This type of harmonic generation is well known in quantum optics [47] and can also be 
realized for elementary excitations in Bose-Einstein condensates [48,49]. As we have shown 
[42], it exists as well for topological coherent modes. 

(8) Parametric Conversion 

When several driving fields are involved, with the frequencies cuj (j = 1,2,...), then 
there exist other possibilities for generating a coherent mode corresponding to the transition 
frequency cj2i- Thus, in the case of two driving fields, with the frequencies uji and UJ2, the 
transition from a mode with energy Ei to that of energy E2 can be activated if UJ1+UJ2 = UI21 
or UJ1—UJ2 — <^2i- In the language of quantum optics [47], this is called parametric conversion, 
up or down, respectively. Generally, for several driving fields, the condition of parametric 
conversion can be written as 

^(±^,) = ^21 . (27) 
j 

An analogous effect is known in quantum optics [47] and for elementary excitations in trapped 
condensates [48,49]. Now, we proved [42] its existence for nonhnear coherent modes. 

(9) Atomic Squeezing 

Atomic squeezing can be defined by means of pseudospin operators, because of which it is 
often termed spin squeezing. For this purpose, one passes from the atomic operators Oj and 
a|, corresponding to the i-th mode and satisfying the commutation relations [oj, a]] = 5ij, to 
the pseudospin operators. In the case of the two-mode condensate, with i = 1,2, one defines 
the collective spin operators 



= a\a2 , S+ = 4ai , ("202 - > 

which satisfy the standard spin commutation relations 

The amount of squeezing can be measured by the squeezing factor 

Q..^. (28) 
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where A'^{Sz) is the dispersion. Atomic operators are connected with the mode amphtudes 
by the equahty < ajaj >= Nc*Cj. Since q = Ci{t) is a function of time, the squeezing factor 
(28) is an oscillating function [12]. 

For the 3- mode condensate, when i = 1,2,3, one introduces the collective angular mo- 
mentum operators 

J_ = \/2 (^a|a2 + 0203^ , J+ = \/2 (alai + 0302) , Jz = 0303 — OiOi , 
with the commutation relations 

[J+, ./_] = 2Jz , [Jz, J±] = ±J± . 
Then squeezing can be characterized by the squeezing factor 

= I < J± > I ' 

similar to factor (28). The generation of atomic squeezing can be useful for spectroscopy. 
(10) Entanglement Production 

To quantify the amount of entanglement produced by a p-particle density matrix pp{t), 
we may employ the measure 

e{p,) = log , (29) 

\\Pp\\-D 

where implies the Hermitian norm over a set T> of disentangled states and means 

the disentangled counterpart of Pp (see details in [50]). For a multimode condensate, we have 

£(pp) = (1 -p)log sup|cjf . (30) 
3 

The mode amplitudes Cj = Cj{t) depend on time. Therefore measure (30) describes the 
evolutional entanglement. The temporal behaviour of the functions Cj(t) can be regulated 

by switching on and off the resonant driving fields. Hence the measure (30) of the produced 
entanglement can also be regulated, which opens the possibility for information processing 
by means of the multimode resonant condensates. 



5 Conclusion 

We have shown that multiple topological coherent modes can be generated in a Bose-Einstein 
condensate, which requires the usage of several driving fields. It is important to emphasize 
that we have accomplished the consideration in two ways, by applying the averaging tech- 
niques and by direct numerical simulations for the Gross-Pitaevskii equation, both ways 
being in good agreement with each other. The multimode condensate presents a novel 
object possessing a variety of interesting properties. From one side it reminds us a multi- 
level resonant atom, so widely considered in quantum optics, but, at the same time, it is 
a multiatomic system, because of which it exhibits imusual nonlinear features. The multi- 
mode condensate displays the following remarkable effects: interference fringes, interference 



9 



current, mode locking, dynamic transition, critical phenomena, chaotic motion, harmonic 
generation, parametric conversion, atomic squeezing, and entanglement production. The 
rich variety of interesting properties, characteristic of the multimode condensates, suggests 
that the latter could find a number of applications. 
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